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The mean field method is applied for analysis of valence electrons in metals. It is shown that at 
low temperatures electrons have two wave-vector distribution patterns. Isotropic distribution refers 
to the first pattern. Anisotropic distribution refers to another pattern, particularly to specific wave- 
vector values occurred nearby the Fermi sphere. It is shown that it is the anisotropy that makes the 
metal obtain its specific superconductor features. 

1. Introduction 

Prior to discussing the matters disclosed in scope of this document, I will give my account of the 
new model Hamiltonian, the Nobel Prize Winner, P. W. Anderson addressed to in his lecture. This 
Hamiltonian is now named, as Anderson Hamiltonian. Hereunder, I'd like to give some extracts from the 
above [1, 2]. 

"To describe the case, I put forward the model Hamiltonian, which is now referred to as Hubbard 
Hamiltonian: 

i j i 

where bij describes the i — to — j node electron jump amplitude and U - repulsion energy of two single- 
centered opposite-spin electrons. Need to say, that no single-centered parallel-spin states are possible. 

It is assumed that there is a local atomic spin denoted by S and existent from the God. It is also 
assumed that this kind of a spin interacts with free electrons on their exchangeable basis. This very 
statement is exactly referred to Kondo Hamiltonian: 



H = ^ £ k n ka + JSS. 



where J is an empirically supported exchange integral, 

k k' ctct' 

- local spin density of impurity conductivity electrons. 

Anderson model is a very simple one to describe the aforesaid electronic mechanism. We introduce 
critically important single-centered U energy and make specific description of the impurity atom by 
applying additional orbital ip d with its occupation number n d and birth operator cj" matched it. Model 
Hamiltonian is formulated as follows: 

H = 2J £k n ka + Un dt n di + E d (n dt + n di ) + ^ V dk c\ a c ka , 

ka ka 

where, apart from free electrons and magnetic term with constant U, we have added tunneling V d k term 
d — k which describes the effect of electron tunneling through the centrifugal barrier to convert orbital 
(p d to one of the Friedel resonances. 
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Before proceeding to other matters of this discussion, I would like to say some words about simplicity 
of the model which is actually more likely superficial than true. Any skills of modeling result in ability 
to give up actual and existent, but not significant aspects of the problem to be of risks for both the 
author and reader. As a matter of fact, the author can give up anything of importance; the reader to be 
in possession of excessively delicate experimental techniques or to be excessively punctual in calculations 
may oversimplify the conceptualized model which is chiefly targeted to demonstration of any specific 
potential. 

As referred to the science style, one of my profound beliefs is that quantum mechanics and statistical 
physics are in principal so simple to make many facts, the Nature present with, look like improbable 
and just the only demonstration of the reasonable interpretable mechanism leaves no doubt in adequacy 
of the interpretation. Moreover, it concerns the cases when any unexpected phenomena are correctly 
predicted; particularly those related to low-density moment life-span correlation, previously described 
orbital-moment quenching, as respects to all the d - level impurities, or negative free-electron exchange 
polarization subject to further discussion. More often than not, such simplified model clears up the 
real nature of the phenomenon, and thereby any 'ab initio' value calculated for the purpose of specific 
cases, which, if they are actually correct, often contain a great number of elements that more likely 
conceal than disclose the truth. Thus, the possibility of providing precise calculation or measurement 
of something may be more likely an obstacle than an advantage, since anything under measurement or 
calculation occurs to be unessential from the standpoint of phenomenon mechanism investigation. After 
all, any ideal calculation just copies Nature, but not explains it". 

How it comes that Hamiltonian is named the model one? Because it is nearly phenomenological - i.e. 
it is not precisely calculable. And if there are a few opportunities to make any calculations, then it is 
fine. This Hamiltonian may be applied for describing the experiment under study. 

2. Anisotropy 

Now, it is time to describe isotropic distribution [3]. We have prescribed function / = /(a), i.e. value 
/ depends on vector a. If value / depends on modulus a of this vector only, the distribution concerned 
is called isotropic, i.e. it may be formulated as / = /(a). This kind of isotropy may be represented 
graphically (see Fig. 1). We will plot a sphere of radius a centered in the origin of coordinates. So, value 
/ will remain equal at any point of this sphere, providing that / = /(a) is the isotropic function. Any 
other / — /(a) function will be referred to the anisotropy one. 
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Fig. 1. Isotropic function. 

Now, we will consider the example of the anisotropic function. We will plot two vectors. One of them 
will be an arbitrary vector a and the other one will be rated as equal, but opposite in its direction — a. 
Two such vectors are shown in Fig. 2. So, if it is appeared that function values fail matching in the points 
concerned, i.e. /(a) 7^ /(— a), this function will be called the anisotropic one. Some exhaustive examples 
of the anisotropic function may be additionally described, but, as a matter of fact, the information 
provided is sufficient for understanding. 
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Fig. 2. Example of anisotropic function. 



3. Fermi Dirac function 

Now, we will consider formulation of Fermy - Dirak function and review its meaning [4] . In the first 
place, we are chiefly speaking about electrons taking account to their large and very large number, for 
example solid-state body electrons. Secondly, electrons are subject to Pauli principle, according to which 
maximum one electron may be in a particular state. This function has one of the following simplest 
formulations: 



where n{e) is an average 
chemical potential. For the function curve, see Fig. 3. 
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Fig. 3. Fermi - Dirac function. 

Fermi - Dirac function shall be derived for the purpose of electron interacting system only. Since 
electrons interact to each other, just an approximate equation may be derived for function n(e). The 
simplest one is mean field approximation in scope of which it is defined that double probability is equal 
to the product of single probabilities. This kind of approximation is known as statistical independence. 
As referred to the mean field approximation, the following equation is obtained: 



n(e) 



e P(e-n) _|_ 1 



(3.2) 



where e is average electron energy which is equal to the sum of kinetic energy e and energy of interaction 
of this electron with other electrons. 

Now, let us remember that in specific representation a particular state of electron is described by 
wave function ^kni^X^iOi where k is a wave vector, n and a - other numbers, which in combination 
with the wave vector make it possible to define the state of electron. The wave vector is linked with an 
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electron momentum by the simple relation p = % k. Now, without going into any details, we formulate 
the equation (3.2) as follows: 

ln l^£k =/3 (e k - M ), (3.3) 
w k 

where w k = n(e k ) i s probability of k vector occupation probability; e k is single electron mean energy; 
e k is electron kinetic energy. By this means as follows: 

£k = £k + X! £ kk'™k'- (3-4) 
k' 

The addend in this formula is the mean energy of interaction of an electron with other electrons. Summand 
s kk > is the energy of interaction of two electrons with wave vectors k and k'. 



4. Model Hamiltonian 

Equation (3.3) contains its nonlinearly unknown function w k . Now, for the purpose of the above 
function, it is necessary to define the electron interaction energy e kk ' . Electrons interact with each other 
under the Coulomb repulsive potential. But there is rather large number of ions and other electrons 
applied for interaction in the solid body. Need to say that it is very hard to calculate the rate of electron 
interaction energy. Therefore, we will use the model Hamiltonian. We will assume that [5] 

£ kk' = 1 <Wk' . 

where c>k+k' is Kronecker delta, / is k and k' = — k wave vector and electron interaction energy. As 
provided in scope of our model, only those valence electrons are repulsive, which are able to surmount 
the crystal at equal opposite direction velocities. 

With the model formula applied (4.1) the average electron energy (3.4) is formulated by the equation 
as follows: 

£ k = £k + ^-k- (4-2) 

According to this formula, the higher k vector electron energy, the more probability w_ k of — k wave 
vector occupation. Hereby, the — k wave vector electron somewhat affects the electron forcing it out of k 
wave vector state. 

We will insert formula (4.2) in equation (3.3). We will gain the following formula: 

In i^L =(3(e k + Iw_ k -v). (4.3) 
w k 

Thus, the equation containing two function values tot and w_ k is produced. 



5. Isotropic and anisotropic electron distribution 

If you ask why it has to do with anisotropy, it may be confirmed by the equation (4.3) which exhibits 
solution of anisotropic function w k subject to condition of w_ k =/= w k . In this equation we will substitute 
vector k for vector — k. If to consider that kinetic energy is the isotropic function, i.e. e_k = £ k ; we will 
formulate the following equation: 

In 1_W - k =/3(e k + Iw k -n). (5.1) 
w_ k 

Equations (4.3) and (5.1) produce the system containing two unknown w k and w_ k . At the same time, 
it is clear enough that probability w k is the composite k vector function, where electron kinetic energy 
is applied as an intervening variable e k : 

w k = w{e k ). (5.2) 

Combined equations (4.3) and (5.1) exhibit their anisotropic solution subject to the condition as 
follows: 

w~ k = w k . (5.3) 
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Using this equation we will eliminate u;_ k from the combined equations (4.3) and (5.1). We will find 
solution of isotropic distribution function by applying the following equation: 

In IZ^JL =/3( £k + / Wk - M ). (5.4) 

There are also some anisotropic distribution functions that fall out of formula (5.3) when specific wave 
vector values are applied: 

w-k + w k . 

Such kind of electron state distribution anisotropy may occur even when no external field is available. 
While applying forms: 

w_k = wi(£ k ), w k = w 2 (e k ), (5.5) 

we may formulate equations (4.3) and (5.1) by the method as follows: 

(5.6) 



In = — (e + w 2 ) , 

W\ T 

In = — (e + wi) , 



w 2 t 

where 



fj, 4 



I ' I 
The following functions remain unknown in the combined equations (5.6): 

Wl=U>l(e) H W2—W2(s). 

If electrons have isotropic wave vector distribution, it is necessary to insert w\ — w 2 = wq in the 
combined equation (5.6). In this case, the equation gained may be formulated by the method as follows: 

e = — in wo . (5.7 

4 w 

This equation states specific dependence of w = w (e) with various temperature values graphically 
represented in Fig. 4 in the form of monotonically decreasing curves. 

If electrons have anisotropic wave vector distribution, probabilities W\ and w 2 in the combined 
equations (5.6) shall be considered as various functions w\ = Wi(e) and w 2 = w 2 (e) subject to energy e. 
To determine these dependences, we will introduce new variables d and s applying the relations as follows 

W2 — w\ = d , Wi + w 2 = 1 + s . (5.8) 

Without loss of generality we will assume that nonnegative difference d of two distribution function values 
W\ and w 2 is d > 0. At the same time, d remains equal to d G [0, 1]. Value s may possess the values within 
the range of — 1 to 1: s 6 [- 1, 1]. We will determine the equalities (5.8), as regards the probabilities w\ 
and w 2 : 

w 1 = ^(l + s-d), w 2 = ^{l + s + d). (5.9) 

We will transform the combined equations (5.6) by applying the formulas (5.9). We will firstly subtract 
specific equation from one of the combined equations and then we will add the equations. As a result, we 
will obtain the following combined equations: 

(1 + rf) 2 -, 2 , d/T , r (1-sf-d* 1 



e 



4d/T > e =Mn , 9 ~" -^(l + s). (5.10) 



{l-df-s 2 ' 8 {l + sf-d 2 2 

The first equation of the combined ones may be easily resolved against s: 



.M-w (i -' ) ' e e-- 1 (i+ "' 
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Fig. 4- Isotropic function of distribution of conductivity electron energy 
at various temperature values t: 1 - t = 0; 2 - r = 0, 25; 3 - t — 0,8. 

As concerns the relations (5.9), probabilities wi and w 2 may be also considered as d functions: wi — 
wi(d), W2 = W2(d). With the second equation of the combined ones (5.10) applied, we may express 
electron energy e in terms of parameter d. Using the dependences produced specific graphs of functions 
wi = tci(f) and w 2 = w 2 (e) may be easily plotted for various temperature values. For the plotted curve, 
see Fig. 5. 
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Fig. 5. Anisotropic function of distribution of conductivity electron energy 
at various temperature values t:1-t = 0;2-t = 0, 25; 3 - r = 0, 8; 4 ~ T = 0, 95. 

The pattern of distribution of electrons by their states depends on their relation between metal 
temperature T and critical temperature: 



At the temperature of T > T c the distribution function w — w{e) is single-valued and satisfies the 
condition (5.3), as respects all the e energy values. At the temperature of T < T c the energy is limited 
by (ei, e 2 ) with function w = w(e) possessing any of three values at every point of the limit, particularly 
wi(e) < u> (e) < w 2 (e)- Being out of the aforesaid limit, the distribution function w = w(e) possesses 
only a single value u>o(e). Thus, equation (5.4) is resolved by applying function = wo( e k) to describe 
isotropic wave vector electron distribution. 
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6. Order parameter 

At T < T c some kind of anisotropic wave vector electron distribution may occur in the narrow layer 
S under Fermi surface e k = fi. This kind of distribution is formulated by 



w k = w 2 (e k ) 



w-k = wi(e k ) 



(6.1) 



subject to e e (ei, e 2 ). Difference d — w 2 — w\ of two anisotropic electron distribution function values 
possesses the largest value d max subject to e = 0, 5. In this case, w = 0, 5 and s = 0. We will determine 
difference d max from temperature r by applying s = in the first equation (5.10): 



2 d max 1 + d max 
= In 



1-dr 



(6.2) 



For dependence curve, see Fig. 6. 




0,5 1 r 

Fig. 6. Electron distribution anisotropy parameter d max , as t temperature function. 



7. Electron T = distribution 

At T = the isotropic distribution function is formulated as follows: 

1 at e k < fi — I , 



1 



(e k — n) at /i — 7<e k <^, 



at £ k > n . 

As for the anisotropic distribution, it is formulated as follows: 

w k = 1 at e k < fi — I , 

w k = 1 , w_ k = or w k = , w_ k = 1 
at fi — I < e k < fi , 
w k = at £k > t 1 ■ 



(7.1) 



(7.2) 
(7.3) 



As provided by the formula (7.2), layer S may be determined under Fermi surface by the inequality (7.3), 
in which the electrons have anisotropic wave vector distribution, i.e. one of the both k and — k wave vector 
states in this layer is free and another one is occupied. For the function curves, see Fig. 7. Apparently, 
electron distribution function obtains its three S layer values. And what is the matter it stands for? The 
answer is in the value of energy, the isotropic or anisotropic distribution electrons exhibit. The electrons 
gain their steady state when they have the lowest energy. 
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Fig. 7. Isotropic and anisotropic distribution of conductivity electrons depending on their kinetic energy 
at temperature t — 0: 1 - isotropic distribution, 2 - anisotropic distribution 



In scope of normalization conditions, the average itinerant electron velocity may be defined by the 
formula as follows: 



where G = 2n s is a number of states in one node and N - mean number of conductivity electrons in 
a crystal. If the distribution function is isotropic, mean electron velocity v gets equal to zero. Formula 
(8.1) may assign specific nonzero electron ordered motion velocity values to some anisotropic distribution 
functions, i.e. these distribution functions are applicable for defining electric current. If there are steady- 
state currents to exist with no external fields available, than such itinerant electron system states shall 
be considered as the superconductive ones [6-8]. 

We will assume that the state of electron gas is described by the anisotropic distribution function 
(6.1) or (7.2). In this case, mean electron ordered motion velocity modulus v may assign any value rated 
from zero to certain v. max . The mean velocity will be equal to zero, providing that free pairs and those 
occupied by wave vectors k and k are chaotically distributed within layer S. If all the states concerned 
are occupied in one half of the layer (this is to say at k x > 0) and free in another half of the layer (at 
k x < 0) , the electrons will gain their maximum ordered motion velocity. The value assigned by the mean 
electron velocity is defined by the nature of initial electron gas state. If the pattern of anisotropic wave 
vector electron distribution is rather steady with respect to small environment variations, the electron 
velocity value will survive for ages. This means that the metal concerned was able to gain its specific 
superconductive characteristics. 

Now, we will calculate the energy the isotropic and anisotropic distribution electrons exhibit. We will 
apply the normalization condition formulated as follows: 



8. Superconductivity. Energy of states 




(8.1) 




(8.2) 



k 



Mean field approximation electron energy takes on the following form: 




(8.3) 



k 
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We will approximate dependence of electron kinetic energy from wave vector k by applying the 
formula as follows: 

£k = -5 , 8.4 

2m 

where m is effective itinerant electron mass. As provided by this formula, any electron kinetic energy 
shall be counted from the band bottom to be also called bottom of conduction band, i.e. e^ = = 0. 

To simplify calculations, instead of k summing we will produce integration by e electron kinetic energy. 
By applying the dependence (8.4) we will obtain the following symbolic equation: 



00 

G]T ... =AN I ... vT 



de . 



where 

GmV2^V fi 2 k 2 _ fr 2 kdk 

/l — o — . s — — ~ , as — . 

2ir 2 h 3 N 2m m 

The upper integration limit may be equal to 00, since the occupational probability of states which energy 
e is specified at the ceiling of conduction band is actually equal to zero. Now, we will formulate the 
normalization condition (8.2) by the method as follows: 

00 

A J w(e)y/7de = 1 . (8.5) 


As for the isotropic distribution electron energy, we will formulate the following formula: 

00 

E (i) = AN J + i / w(e)j w(e) ^JT de . (8.6) 


If the isotropic electron distribution function applied at T — is formulated according to (7.1), the 
equations (8.5) and (8.6) take on the following form: 

n-i n 



A J yfe de + A J ii_£ ^/J de 





H-I — n 

=AN J ^+I/j v Td e + 4f J^ 2 -e 2 )^7de. 
n-I 
Since we apply the small parameter 

A = — . 
e F 

where 

2/3 

e F 



2 A 



- refers to Fermi energy, we will define that the chemical potential and isotropic wave vector distribution 
electron energy at T = take on the forms as follows: 
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We will assume that the anisotropic wave vector electron distribution at T = is defined by the 
function as follows: 

1 at £k < fi — I , 



1 at a — I < eu < u, , k x > 

(8.8) 

at jj, — I < s k < jj, , k x < 
at e k > \i . 

As provided by the above formula, only one half of the k x > states may be referred to as the occupied 
ones to occur in layer 5* above the Fermi surface, which thickness 5k is proportional to interaction 
parameter I. For anisotropic distribution pattern, see Fig. 8. 




Fig. 8. Anisotropic distribution function at T = 0. 
Here, the normalization condition gives rise to the following equation: 

n-i 



A J yfe de+^A j yfe de = 1 . 



[i-I 

As for electron energy, it may be calculated by the formula as follows: 



J (e + 7T l) yfe de + ^ A T / : K •""./_-. 



E Q = AN 

o n-I 
As provided by the above calculation, the following formulation is obtained: 



(8.9) 



If occupied and free state pairs that match specific wave vectors k and k will be distributed within 
layer S by any other way, the chemical potential and electron energy rating will remain the same. 

The difference of electron energy values (8.7) and (8.9) will be formulated by the equation as follows: 



NI 2 

8e F 



> 0. 



Thus, we get to the conclusion that the state of itinerant electrons described by the anisotropic distribution 
function is the primary one - i.e. the electron system specified in this condition is of the lowest energy. 
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Considering for the aforesaid about the anisotropic electron energy distribution we will plot the pattern 
of superconductive state, as shown in Fig. 9. 



w 








e — fi 



Fig. 9. Anisotropic distribution of energy e conductivity electrons 
subject to the lowest energy E at the temperature of t — 0.8. 



9. Maximum superconductivity electron velocity at T = 

Now, we will find ordered electron motion velocity in the state described by the distribution function 
(7.2) at the temperature of T = 0. For this purpose, we well substitute the wave vector sum specified in 
the formula (8.1) for the following integral: 



— =■ / kwird 6 k . 

',7T 3 mN J 



(9.1) 



If T = 0, the inner and outer radii of layer S shall be respectively equal to as follows: 

1 



fei = — y/ 2m(/j-/) , 



fc 2 = — J 2 m a . 
h 



In this case, the layer thickness is defined by the following formula: 

I 



5k = fc 2 — ki — 



m 
2e F 



(9.2) 



If the states are filled with electrons in one half of layer S, but another one remains free, the rate of 
ordered electron motion velocity exhibits its maximum value 



3/ 

4 y/2meF 



(9.3) 



10. Electric current-forced superconducting state destroy effect 

Superconducting state of the itinerant electrons is destroyed in the events when the current running 
over metal items exceeds its particular critical value (Silsbee effect). We will assume that specific homogeneous 
electric field with tension E is produced inside the metal. Under the effect of this field the itinerant 
electrons will execute their ordered motion at the average velocity u which direction agrees with another 
one that affecting force electron — eE. In this case, electron state distribution function may be defined 
by the equation (3.3), where electron energy depends on wave vector k in the following manner: 

= e k — h k u + I w -k ■ (10-1) 
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It is no easy matter to make exact solution of the equation (3.3). Therefore, we will consider its approximate 
solution only. We will assume that distribution function has the form as follows: 



»k = /(k- k Q ) , 

where = /(k) is the solution of the equation (3.3), providing that u = 0; 



(10.2) 



k 



m u 

% 



Function to^ = /(k) is equal to a unity almost at all the points occurred inside Fermi sphere: k < kp, 
excluding the points at surface S of that sphere. If occurred outside Fermi sphere, function = /(k) 
almost everywhere is equal to zero. Field D containing nonzero function (10.2) is limited by the sphere 
with radius kF, which center C is displaced off the origin of coordinates O by vector k . In other words, 
the wave vectors that agree with occupied electronic states belong to field D. The theory under discussion 
is applied to superconductivity to be due to interaction of electrons, which wave vectors k and — k belong 
to spherical layer S with its radius equal to kp and thickness - to 5k. No anisotropy is created with the 
electrons distributed over their wave vectors, when the rate of displacement k of field D is so large that 
one half of layer S finds itself beyond this field. For arrangement of field D and layer S, see Fig. 10. 




Fig. 10. Displacement of Fermi sphere under influence of electric field. 

Point A belongs simultaneously to the displaced Fermi sphere and to the inner surface of layer S. 
Therefore, 

AC = k F , AO = k F ~5k. 
We will apply the rule of Pythagoras for AOC triangle. Now, we will gain the following formula: 

k F = k 2 Q + (k F - 5k) 2 . 

It appears from this equation that superconductivity-force displacement k will be formulated as follows: 

k Q — \J 2kp 5k . 

This formula may be transformed to: 



kn — 



V 2m& 



where 5e is energy width of the layer S. In this connection, the average current velocity will be calculated 
by the following equation: 



2 5e 



(10.3) 



By this means that the superconductive state of conduction electrons will be destroyed when the external 
electric field makes them moving to the same direction and produces the current, which density exceeds 
the following value: 

j K p = enu. 
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We will specify the relation of maximum superconductive current velocity v max to the least current 
speed u destroying the superconductive state: 

^ = -,/^. (10.4) 
u 8 V e F x 1 

It is clear enough that the above relation is much less than unity v max -C u. 

On cutting off the electric field and after thermodynamic equilibration the electrons shall transform to 
their superconductive state distributing as provided by the formula (7.2). Thereafter, the average electron 
velocity shall drop down to value v max and superconductive current of maximum density j max will run 
over the metal. 



11. Mean energy dependence of kinetic energy 

The rate of mean electron energy dependence of its kinetic energy is defined by formula (4.2). 
As may be inferred from the above formula, the electron energy with wave vector k depends on whether 
the — k wave vector state is free or occupied. Electron energy e may be specified by the e kinetic energy 
functional form as follows: 

e(e) = e + 1 wi(e) . 
or the pattern of this function at various temperatures, see Fig. 11. 




Fig. 11. Mean electron e energy dependence of its e kinetic energy 
at various temperature values t: 1 - r = 0, 2 - t — 0,5, 3 - r = 0, 8. 



12. Conclusion 

Thus, in scope of the theory discussed the microscopic superconductivity is due to anisotropic wave 
vector electron distribution. Normally, anisotropy is caused by electron repulsion effect in the k and — k 
wave vectors states. If applied to the matter of more complex nature, model Hamiltonian is formulated 
as follows: 

e kk' = I s k+k' - J ^k-k' i 

where J is attractive energy of two electrons in the conditions equal to wave vectors k = k'. This 
Hamiltonian is calculated in [9]. 
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